Final spin of a coalescing black-hole binary: an Effective-One-Body approach 
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We update the analytical estimate of the final spin of a coalescing black-hole binary derived 
within the Effective-One-Body (EOB) approach. We consider unequal-mass non-spinning black- 
hole binaries, ft is found that a more complete account of relevant physical effects (higher post- 
Newtonian accuracy, ringdown losses) allows the analytical EOB estimate to "converge towards" the 
recently obtained numerical results within 2%. This agreement illustrates the ability of the EOB 
approach to capture the essential physics of coalescing black-hole binaries. Our analytical approach 
allows one to estimate the final spin of the black hole formed by coalescing binaries in a mass range 
( = mim2 / {mi + ni'if' < 0.16) which is not presently covered by numerical simulations. 

PACS numbers: 04.25.Nx, 04.30.-w, 04.30.Db 
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I. INTRODUCTION 



The dimensionless spin parameter a — a/M — J/M^ 
of the black hole formed by the coalescence of a black-hole 
binary is a complicated function of the initial masses and 
spins of the two constituent black holes. [We consider 
inspiralling systems circularized by radiation reaction]. 
The function a(mi, TO2, Si, S2) is a useful diagnostic for 
comparing analytical approaches to the dynamics of coa- 
lescing binaries to the results of three-dimensional numer- 
ical simulations. Here we shall consider the simple case 
of non-spinning binaries, where d becomes simply a func- 
tion of the symmetric mass ratio v — 77117712/ {mi + 7712)^ 
(which satisfies < < 1/4) . 



As far as we know, the first estimate of a was made 
in 2000 [T] , on the basis of a new analytical approach to 
the general relativistic two-body dynamics, the Effective- 
One-Body (EOB) approach 0, at a time where no reli- 
able numerical simulations of coalescing black hole bina- 
ries were yet available. This estimate was a = 0.795 for 

= 1/4 (i.e. for the equal mass case), and was based 
on estimating the ratio J/M^ at a "matching radius" 
''match — 2.85M where the two-black-hole system was re- 
placed by a unique, ringing black hole. This estimate 
used a 2.5 post-Newtonian (PN) accurate description of 
the dynamics down to rmatch, and neglected the energy 
and angular momentum losses during the ringdown. In 

2001, a combination of full numerical simulations with 
a "close-limit" approximation Q describing the ringing 
final black hole led to a similar estimate, namely a ~ 0.8 
for 1/ = 1/4 [J]. This estimate was revised downwards in 

2002, to a ~ 0.7 @, when a better {Lazarus-ha.sed) way 
of bridging the far- and close-limit approaches indicated 
that the angular momentum losses were larger than pre- 
viously estimated. Indeed, Ref. Q had estimated the to- 
tal angular momentum loss after crossing the Last Stable 



Orbit (LSO) ^ to be around 2%, and the corresponding 
energy loss to be around 3%. [The EOB approach had 
estimated an energy loss beyond the LSO of about 1.4%, 
half of it emitted during the plunge and the other half 
during the ringdown [^]. By contrast, Ref. 0] estimated 
the angular momentum loss to be around 12%. More re- 
cently, an update of the analytical EOB estimate JJ using 
the now available 3.5PN accurate description of the two- 
body dynamics showed that increasing the PN accuracy 
of the dynamics (from 2.5PN to 3.5PN) had the effect 
of decreasing the final spin parameter from d = 0.795 to 
d — 0.77 (again for u = 1/4, and again when neglecting 
the angular momentum loss below the matching radius 

'"match) • 

The estimates just recalled belong to a "prehistoric" 
era where, for a variety of reasons, numerical simu- 
lations did not exhibit a very convincing convergence 
among themselves, nor towards analytical results (see, 
however, 5] and 8]). This era has recently ended thanks 
to remarkable breakthroughs in numerical relativity. Dif- 
ferent groups have finally succeded in numerically simu- 
lating the merger of two black holes of comparable masses 
mi and 777,2, possibly with spin, and their results exhibit 
convincing internal convergence, and a nice mutual con- 
sistency [a, [13, [nl, [11, [11, til [11 [B [IS [ii, [ii, [13 . in 

the particular case of non-spinning black holes with equal 
masses, ttii = m2, different groups now agree on the value 
d ~ 0.69 (which is within the span of the "prehistoric" 
estimates recalled above) . The most extensive analysis of 
the final angular momentum of coalescing black holes (in 
the non-spinning case) to date has been carried out by 
Gonzalez et al. [1^1 : they have considered a large sample 
of unequal-mass systems, corresponding to a symmetric 
mass ratio v varying in the range 0.1613 < v < 1/4, 
and have accurately determined the variation of d with 
v within this range. The aim of the present paper is 



^ Also known as the innermost stable eircular orbit (ISCO). 
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to generalize and update the analytical EOB estimates 
P, [3] of a(^) recalled above, both by explicitly consid- 
ering general values of v and by improving the previous 
EOB treatments of the physical ingredients which are 
crucial in determining the value of a. Our motivation for 
this study is two-fold: on the one hand, we wish to see to 
what extent numerical results can be reproduced by ana- 
lytical (EOB) estimates, and on the other hand we wish 
to understand, on this example, what are the physical 
ingredients which are crucial in accurately determining 
the plunge, merger and ringdown dynamics of coalescing 
black holes. 

Let us recall that the EOB approach to the general 
relativistic two-body dynamics is a non-perturbatively re- 
summed analytic technique which has been developed in 
Refs. [H, i, 0, iH m [23(. This technique uses, as ba- 
sic input, the results of PN theory, such as: (i) PN- 
expanded equations of motion for two point-like bodies, 
(ii) PN-expanded radiative multipole moments, and (iii) 
PN-expanded energy and angular momentum fluxes at 
infinity. For the moment, the most accurate such re- 
sults are the 3PN conservative dynamics [l^, , and the 
3.5PN energy flux [H, [13, ill . Then the EOB approach 
"packages" this PN-expanded information in special re- 
summed forms which extend the validity of the PN results 
beyond the expected weak-field-slow- velocity regime into 
(part of) the strong-fleld-fast-motion regime. The aim 
being to use the EOB approach for analytically describ- 
ing the last inspiralling orbits, the transition from inspi- 
ral to plunge, and the plunge itself down to a "matching 
radius" rmatch small enough to allow one to match there 
the plunge waveform to a ringdown one. 

The basic new ingredient used below to improve the 
previous EOB estimates of d is an approximate treat- 
ment of the energy and angular momentum losses dur- 
ing ringdown. These losses were neglected in [H, 0|- 
Our approximation will consist in estimating these losses 
by rescaling (proportionally to v'^) the losses obtained 
by numerically studying the test-mass limit a la Regge- 
Wheeler-Zerflh [29i, 30J. We shaU also study the effect of 
the radiation reaction force during the plunge. Instead 
of using (as in [l|) a naive analytic continuation of the 
radiation reaction appropriate to the inspiral phase, we 
shall also use another Pade resummed radiation reaction 
force, which does not assume the validity of Kepler's law 
j-j2^3 ^const. during the plunge [1^. Indeed, the recent 
Regge-Wheeler-Zerilli-like study of the waveform emit- 
ted by plunging test-masses [23, [30] has shown that this 
modified radiation reaction stays closer to the "exact" 
gravitational wave angular momentum flux computed a 
la Regge-Wheeler-Zerilli. 

The paper is organized as follows. In Sec. [IT] we review 
the non-perturbative construction of the two-body dy- 
namics incorporating radiation reaction effects while we 
devote Sec. IIIII to the presentation of our results. Some 
conclusions are presented in Sec. IIVI We use geometric 
units G — c = 1. 



II. EQUATIONS 

In this section we recall the non-perturbative construc- 
tion of the two-body dynamics including a radiation re- 
action force. We take advantage of the most complete 
PN results, i.e. we work at 3PN for the conservative 
part of the dynamics and at 3.5PN for the radiation 
damping. In the EOB framework, the complicated PN- 
expanded relative dynamics (in the center of mass frame) 
of the binary system of masses mi and m2 is mapped^ 
into the simpler geodesic dynamics of a particle of mass 
fi = TOim2/(mi 4- 7712) moving in some effective back- 
ground geometry (in Schwarzschild gauge) 



ds^ 



-A{r)dt^ + B{r)dr^ + {d9^ + sin^ Odip"^) . (1) 



Here, and below, we work with the dimensionless reduced 
variables r = R/M and t = T/M, with M ^ nii + 
7712] (r, 0, 1^) are polar coordinates in the effective problem 
that describe the relative motion.The coefficients of the 
effective metric at the 3PN approximation |2lj read 



A^^^(r) = l-- + 



2iy 



94 41 . 
y"32"^ 
61' 



v 



(2) 



{BAY'-'' (r) = D-''-'' (r) = 1 - ^ + 2{2,v - 26)- . (3) 

Note that we work with the PN expansion of the quan- 
tity D{r) = B{r)A(r)^ rather than with the g^r metric 
coefficient B{r) = D(r)/A{r). [Recall that B{r) is equal 
to 1/A{r) in the Schwarzschild case (which corresponds 
to the test- mass limit ^ 0).] 

Though the EOB packaging of the complicated orig- 
inal PN dynamics into the much simpler metric coeffi- 
cients A{r),D{r) already represents an efficient resum- 
mation of PN-expanded results, it is quite useful, espe- 
cially at 3PN, to further resum the 'packages' v4(r), D{r). 
Indeed, a general feature of the EOB philosophy is to 
smoothly connect the EOB structures to their v 
limit [21, '2^. Here, this can be done by replacing the 
(Taylor expanded) metric coefficients given by Eqs. (|2|)- 
([5]) with suitable Pade approximants. The simplest, and 
most robust, choices consist in using, as metric func- 
tions, the following definitions A{t) = Pg^]^'^^'^] and 
D{r) = Pg [D'^^'^] ^. These Pade approximants are used 
to ensure the following two facts: (i) the function A has a 
simple zero for a positive value of r (like ^schw = 1 ^ 2/r) 
and (ii) the function D stays positive while r decreases 
(like Dschw = !)• The Pade resummation of A is useful 
for ensuring the existence and i^-continuity of a last sta- 
ble orbit (LSO), as well as the existence and ^/-continuity 
of a last unstable orbit, i.e. of a i/-deformed analog of a 



^ modulo some quartic-in-momenta additional terms; see below. 
3 We recall that Pg^ID^PN] = 1/(1 + gj^^a _ 2(31, _ 26)vu^) where 

u = 1/r , while P3 [A^^^^] is explicitly given in Eq. (4.40c) of 

Ref. [2I. 
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light ring (LR) . We recall that the LR corresponds to the 
circular orbit of a massless particle, or of an extremely 
relativistic massive particle, and is technically defined 
by looking for the maximum of A(r)/r'^ , i.e. by solving 
{d/dr) {A{r) jr^ ) = 0. The Pade resummation of D is use- 
ful to ensure that the orbital frequency — dipjdt has a 
clear maximum at (approximately) the EOB j/-deformed 
light ring and then drops to zero (cf Fig. [1] below for 
i.= l/4). 

In the EOB approach one splits the general relativistic 
relative dynamics of a binary system into a conservative 
part, determined by the EOB Hamiltonian defined below, 
and a non-conservative part related to the loss of angular 
momentum through gravitational radiation. The EOB 
Hamiltonian (divided by /i) is given by 



tonian 



(4) 



where /fcff denotes the so-called "effective Hamiltonian" 
(describing the geodesic dynamics of the "effective" test- 
mass /i), originally written as pl| 



(5) 



Here zg = 2v{^ - iv), H = H/^i, H^s = H^s/tJ^. 

= P^/{fiM), r = R/M and Pr denotes the conjugate 
momentum to r. 

Following what we did in the test-mass limit case (29l . 
Isoj , the relative dynamics is somewhat more conveniently 
described by replacing the Schwarzschild-like radial vari- 
able r by the EOB generalization of the Regge- Wheeler 
tortoise coordinate , defined by integrating 



dr 



1/2 



(6) 



where we recall that B — D/A. One then needs to re- 
place the r-conjugate momentum pr by the r*-conjugate 

momentum p^. , such that Pr,dr^ = Prdr, i.e. Pr, = 
1/2 

{A/By'pr. The reason for using this transformation is 
that Pr, has a finite limit when r tends to the zero of 
A{r) ( "j/-deformed effective horizon"), while pr diverges 
there. Actually, as we stop evolving the dynamics around 
the i/-deformed light ring (i.e. before reaching the zero 
of A{r)), this change of variables is not really necessary. 
It is, however, convenient because it magnifies the ra- 
dial axis in a crucial region, and prevents any excessive 
growth of the radial momentum during the plunge. [Let 
us mention in passing that the same kind of coordinate 
has been used in [31j] for studying gravitational pertur- 
bations of non-rotating relativistic stars.] 

Neglecting (as it is consistent in a 3PN correction term) 
the square of the factor B/A entering the z^pf/r'^ term, 
this leads to the following form for the effective Hamil- 



\ 



Pr, 



(7) 



Hamilton's equations for {r,ip,pr,,Pip) then read 
Ap^ 



dip 
'dt 




(8) 

(9) 
(10) 



A' 2A 



(11) 



where A' = dA/dr. In these equations the extra term 
represents the non conservative part of the dynamics, 
namely the radiation reaction force. 

During the quasi-circular inspiral, a rather accurate ex- 
pression for JF^ is the following Pade-resummed form [3^ 



32 7/3 /Dis(^n;t^) 



(12) 



which is expressed in terms of the PN ordering parameter 
Vfi = fl-^/^. In this expression, the function /dis denotes 
the "factored flux function" of Ref. [Ill, scaled to the 
Newtonian (quadrupole) flux [see Eq. (4.6)-(4.8) there]. 

Ref. (1|] assumed that the analytical continuation of the 
expression might still be a sufficiently accurate de- 
scription of radiation reaction effects durin g; th e plunge. 
On the other hand, the authors of Ref. [23| pointed 
out that the expression (fT^ assumed the continued va- 
lidity of the usual Kepler law ^ U^r^ = 1 during the 
plunge. [This is why we label the expression p2p with a 
superscript K, for Kepler.] They, however, emphasized 
that the Kepler combination K = fi^r"^ significantly de- 
viates from one after the crossing of the LSO, to become 
of order of 0.5 at the (effective) light ring. Ref. [2^ went 
on to argue for a different expression for the radiation re- 
action, say (without any superscript), that does not 
assume Kepler's law. This new expression reads 



l-v /v^^^ 



(13) 



^ modulo a factor 4> taken into account below. 
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FIG. 1: Time evolution of the orbital frequency Q for 1/ = 
1/4: the maximum occurs near the z/-deformed light ring at 
rLR(l/4) « 2.316. 



TABLE I: A sample of the numerical data of Fig. [2] From 
left to right the columns report: the symmetric mass ratio v, 
the final dimensionless angular momentum aJ^^' from Ref . [l^ , 
and our best estimates (with 3PN+3.5PN dynamics) with 
{a^^,M§^/M) and with (aBH,MBH/M) . 
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where v^p = ^r^^ and r^^ = rib^ where the function ■)/) is 
defined as in Eq. (22) of Ref. [23|^. Note that the essential 
difference between the two possible expressions for the ra- 
diation reaction is that cx fi^/^, while !F^p cx Sl^r*. 
See Ref. jMi] (notably Fig. 2 there) for a detailed com- 
parison of these two analytical representations of radi- 
ation reaction to "exact" numerical results during the 
plunge. In both possible expressions for the radiation re- 
action, our current "best estimate" of / is obtained by 
Fade approximating the currently most complete post- 
Newtonian results, namely the 3.5FN ones [26, 27, 28]. 

In the forthcoming analysis we shall compare (and con- 
trast) the relative dynamics and the related final black 
hole spin using various FN-accuracies for the EOB dy- 
namics (2FN Hamiltonian -I- 2.5FN radiation reaction, 
versus 3FN Hamiltonian -I- 3.5FN radiation reaction), as 
well as the two different expressions for radiation reaction 
briefly discussed above. 



III. RESULTS 

The computation of the relative dynamics needs two 
separate steps: (i) to initialize the system ([5])- pT|) and 
(ii) to integrate it in time. The initial condition for the 
relative dynamics is given in a standard way, notably by 



^ The quantity is introduced to simplify some expressions, be- 
cause it is such that the Kcpler-looking law 0?r^ = 1 holds, 
without correcting factor, during the inspiral (i.e. above the 
LSO). 



specifying a non-zero initial value for p,.^ according to 
the FN order that is being used. Our implementation 
follows Eqs. (4.16)-(4.21) of Ref. 1] and Eqs. (4.8) and 
(4.10) of Ref. 0] [see also Eqs. (9j-(13) of Ref. [29|] and 
does not need to be discussed explicitly here. Let us only 
mention that, for the initial relative separation that we 
shall take, namely ro = 10, the leading post-adiabatic 
approximation is sufficient for getting a smooth quasi- 
circular inspiral (without noticeable eccentricity). 

The orbital frequency O develops a maximum at ap- 
proximately the location of the last unstable EOB cir- 
cular orbit defined by the condition (v4(r)/r^)' = 0. As 
already mentioned above, in the test-mass limit (z^ <C 1), 
this condition defines the light ring r — i. When v ^ Q, 
we shall refer to the solution of (^(r)/?'^)' = as the 
i^-deformed light ring (LR): ri^^{v). 

In the 1/ <C 1 limit, it was realized long ago {3^ that 
the crossing of the light ring by a test particle corre- 
sponds to triggering the black hole quasi-normal modes. 
For related reasons (discussed in [s^l), in the comparable- 
mass case, the crossing of the i/-deformcd light ring cor- 
responds to an abrupt change of description: before this 
crossing one can still describe the two black holes as two 
point masses with EOB relative dynamics, while after 
this crossing one can replace the binary system by a sin- 
gle distorted black hole (as in the "close-limit" approxi- 
mation of colliding black holes Q). In other words, the 
EOB approach estimates the full waveform by match- 
ing at Tniatch — ''LR the inspiral 4- plunge waveform 
computed from the EOB dynamics to a superposition 
of quasi-normal modes (QNMs) describing the ringdown 
of the final distorted black hole. 

Refs. [l|, 0] then estimated the mass and angular mo- 
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FIG. 2: Dimensionless spin parameter of the final blacli hole: 
convergence of Effective-One-Body estimates towards Numer- 
ical Relativity ones as one increases the post-Newtonian ac- 
curacy and takes into account the ringdown losses. See text 
for discussion. 



mentum of the final black hole by the (EOB) energy 
and angular momentum of the binary system at the 
matching point^ ?'match- This procedure, with the choice 
''match = fhRii^), givcs the foUowing estimate of the mass 
and spin parameter of the final black hole 



pLR 



(14) 



This leading- order estimate of obhC'^) is plotted for dif- 
ferent post-Ncviftonian approximations in Fig.[2]and com- 
pared with the numerical relativity data of Ref. (l6j 
(black solid line). The labelling of these leading-order 
"light ring" analytical estimates indicates both the choice 
of the PN accuracy for the dynamics ("2PN-LR" or 
"3PN-LR"), and the choice of a specific radiation reac- 
tion expression, namely Eq. ('Rad Reac oc il^/^') or 
Eq. US]) ('Rad Reac oc n^r^'). 

It is evident from Fig. [2] that these "LR" approxima- 
tions overestimate the "actual" result. It is also clear 
that using a 3PN-accurate dynamics, instead of a 2PN 
one, goes in the good direction. 

As we were mentionning in the introduction, the ba- 
sic physical ingredient which is lacking in these simple- 
minded LR estimates is the loss of angular momentum 
during the ringdown phase (the importance of this loss 
was first emphasized in Q). Here we shall use, as 



next approximation beyond the leading-order estimates 
(|14p . an approximate treatment of the ringdown losses 
of angular momentum (and energy) which takes advan- 
tage of the corresponding results in the small mass case. 
Recently, by combining EOB and Regge-Wheeler-Zerilli 
techniques, Ref. [1^ computed the gravitational wave 
signal (decomposed in multipoles) from the late inspi- 
ral, plunge and ringdown in the small v limit. In this 
calculation the dynamics of the source {r{t),ip{t)) is de- 
termined together with the waveform ^^^(u), and this 
allows one to relate the retarded time u, in terms of which 
the waveform is expressed, to the dynamical time t, en- 
tering the equations of motion ([8|)- pT|) . In particular, 
this allows one to define precisely the retarded time, say 
UhR, corresponding to the crossing of the light ring, and 
thereby the part of the waveform which corresponds to 
the ringdown (in the EOB sense). It was also verified in 
Ref. [29| that the ratio 'i'lmiuj/fj, had a universal limit 
as I' — > 0. By integrating over u > mlr quadratic expres- 
sions in ^(m{u)/n and its derivative [29|, and summing 
over the various multipoles up to ^ = 4, one obtains the 
following numbers for the energy and angular momentum 
losses during the ringdown, in the limiting case v 1: 

(M//L(2)i?""Sdown ^ 0.2448 and ^-2 jringdown ^ ]^ gggg ^ 

Though we know (and discuss below) that the propor- 
tionality of these losses to fj? oc v'^ is only strictly valid 
when <^ 1, our proposal here is to define a next-to- 
leading order approximation (beyond the leading-order 
one (|14p ) by simply rescaling those losses (proportionally 
to ly'^y up to any finite value of ly. In other words, let us 
consider that the angular momentum and energy losses 
due to the ringdown phase {u > ulr ) are approximately 
given by 



E- 



ringdown 
scaled 



J. 



ringdown 
scaled 



(i/) ~ 0.2448 J/^M 
(i/) ~ 1.3890 i/^M^ 



(15) 
(16) 



so that we can define the following next-to-leading order 
approximations to the mass and spin of the final Kerr 
black hole: 



MBH(t^) 
-^BH(t^) 



calcd 
ringdown / \ 



(17) 

(18) 
(19) 



We plot in Fig. [5] the next-to-leading order estimate 
(fT9)) for both choices of radiation reaction and for dif- 
ferent post-Newtonian approximations (they are labelled 
by, e.g., '2PN: Rad Reac oc 9^1'^'). It is clear on Fig. [2] 
that these next-to-leading order estimates are closer to 



Actually, Ref. 01 used as end point of the EOB evolution a point 
determined by the breakdown of certain adiabatic conditions. 



^ We use here the fact that such simple-minded v rescalin gs h ave 
been found to be (surprisingly) rather accurate, see e.g. j34| for 
the energy loss in the head-on collision of two black holes, and 
its comparison to the test-mass limit |33|| . 
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the numerical relativity ones. It is also clear that the 
use of the 3PN-accurate dynamics improves significantly 
the results, compared to the 2PN-accurate case. Actu- 
ally, we see that the 3PN-I-3.5PN order next-to-leading 
estimates (IT9)) closely "bracket" the numerical relativity 
result. Using the (originally proposed) "Kepler- type" ra- 
diation reaction (|I2p leads to a spin parameter which is 
slightly smaller than the numerical relativity one, while 
using the more recently proposed radiation reaction (jI3[) 
leads to a spin parameter which is slightly larger than the 
numerical one. Note also that our analytical predictions 
are computed for all values of v within the full range 
< < 0.25. This allows us to predict (e.g. by using 
the 3PN: Rad Reac oc Vi'/'^ curve, or some combination 
of the two 3PN curves) the angular momentum of the fi- 
nal black hole for the range of values v < 0.16, which has 
not been explored (yet) by numerical simulations ( but in 
which the analytical approximation should be rather re- 
liable). For the sake of comparison, we list in Table U our 
best (bracketing) 3PN-I-3.5PN numbers for the dimen- 
sionless spin parameter, together with a selected sample 
of the numerical relativity data. We mention, in pass- 
ing, that had we estimated the spin parameter as the 
ratio (P^^ - ^scafod°'^'')/-^LR (i-e- neglecting the energy 
loss due to ringdown) , and worked at the 3PN level with 
our (a priori preferred) radiation reaction T^p, we would 
have obtained values even closer to the numerical relativ- 
ity ones: e.g. a(0.25) ~ 0.6804, a(0.1825) ~ 0.5281,and 
a(0.1613) ~ 0.4757. This (probably partly accidental) 
agreement illustrates the fact that our current (approxi- 
mate) analytical framework has already captured most of 
the correct physics, and that small variations in the im- 
plementation of the EOB approach can probably lead to 
an excellent agreement with numerical relativity results. 



IV. CONCLUSIONS 

In conclusion, we have improved previous analytical 
estimates, derived within the Effective One Body (EOB) 
framework, of the final spin of an unequal-mass coalesc- 
ing black-hole binary. Our improvements consist in tak- 
ing a more complete account of the most relevant physics: 
we have used higher PN accuracy, we considered sev- 
eral ways of modelling the radiation reaction during the 
plunge, and, most crucially, we took into account the an- 
gular momentum and energy lost to gravitational radia- 
tion during the ringdown. Our final results differ by less 
than 2% from the recent numerical relativity estimates 
of [l^ (see Table IJ). This nice agreement shows, in 
our opinion, the ability of the EOB approach to capture, 
qualitatively and quantitatively, the essential physics of 
the plunge and merger of black-hole binaries. 

We did not try here to further reduce the small re- 
maining difference between analytical and numerical re- 



sults^. Our goal here was mainly to exhibit, in a simple 
case study, how the inclusion of more and more phys- 
ical effects in the EOB approach led to a nice, mono- 
tonic convergence towards a numerical relativity result. 
In separate investigations [3Q . i35] we shall illustrate how 
the EOB framework can also nicely converge towards the 
gravitational waveform. The aim of these studies is to 
understand which parts of the physics included in the 
EOB method must be more precisely modelled to yield 
accurate representations of the various physical observ- 
ables of merging binary black holes. Indeed, the general 
philosophy of the EOB approach is that this (resummed 
perturbative) analytical framework contains several flexi- 
bility parameters which can be determined by fitting EOB 
predictions to some non-perturbative data, such as nu- 
merical relativity simulations, or, possibly, actual obser- 
vational data. An example of a flexibility parameter is 
the coefficient 05(1^) parametrizing presently uncalculable 
4PN (or higher) additional contributions +a5{i^)/r^ + . . . 
to the crucial "radial potential" A{r) in the effective met- 
ric Eq.([T|). Ref. Q exemplified how the parameter 05 
(such that a5{v) = a^v) could be fitted to numerical data 
(for initial configurations). 

We are aware of the rather coarse nature of the approx- 
imation used above for estimating the ringdown losses. 
In fact, the EOB approach itself provides a better, and 
more consistent, way of estimating these losses. Indeed, 
by matching, at rmatch — '"lr, a- post-Newtonian im- 
proved plunge waveform (£ — 2, m — ±2) (from the 
EOB 3PN-I-3.5PN dynamics) to a superposition of quasi- 
normal modes (QNMs) of a Kerr black hole of mass and 
angular momentum given by the energy and angular mo- 
mentum of the relative dynamics at the ;/-deformed EOB 
light ring, we can analytically determine the amplitude 
of the ringdown waveform. Then, from this waveform we 
can analytically estimate, within the EOB approach, the 
ringdown losses. In view of the many delicate issues con- 
nected with this matching procedure (see [13, [H, \3T\). 
we leave to future work a detailed discussion. Let us, 
however, quote some preliminary results that we have 
obtained. In the = 1/4 case, and considering, for sim- 
plicity, only the contribution of the quadrupole {£ — 2, 
m = ±2) part of the waveform, this matching proce- 
dure gives JEOTmrched,2,±2(l/4) - 0.0899 for the angu- 
lar momentum carried away by^ = 2,r7i = ±2 gravi- 
tational waves (GW) during ringdown. This happens to 
be in good agreement (relative difference < 4%) with our 
naively scaled estimate: J™fg'^°™(l/4) w 0.0868. How- 
ever, this good agreement is partly accidental. Several 
complicated effects go here in various (probably com- 
pensating) directions: (i) the matching can under- or 



Nor do we wish to conclude from Table I that is a more 
accurate expression for the radiation reaction than JP^. One 
needs to consider the effect of higher PN contributions before 
reaching any conclusion. 
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over-estimate the amphtude of the ringdown signal, (ii) 
the exact scaling with v is not exactly proportional to 

^ (iii) in the = 1/4 case, the GW signal is dom- 
inated by the quadrupole £ = 2, m = ±2, (iv) for 
smaller values of v (when the reflection symmetry is 
lost) the higher multipoles provide significant contribu- 
tions, and, indeed, our scaling estimate was based on 
the small v limit and included all the multipoles up to 

= 4. [As an example of the importance of higher 
multipoles for smaller i/'s let us mention that we indeed 
found, for = 0.1 a quadrupole-only matched value of 
•^ESBrt'ched,2.±2(0-l) ^ 0-00631 which is roughly half the 

scaled value J™ld°™(0-1) ~ 0.0139]. 

As a further remark, we point out that our analyti- 
cal results above suggest that a simple quadratic expres- 
sion h'nx.iy) = CL\v + a2v'^ should provide a reasonable fit 
to the data. And indeed, fitting the data of [l^ with 
this function, provides a good fit^ when a\ — 3.27690 
and a-i = —2.11405. Moreover, this fit is found to stay 
close to our best analytical predictions for the range 
V < 0.16 not covered by numerical simulations. For 
instance. For i/ = 0.1, this gives afit(O.l) ~ 0.3066, 
to be compared with, say, a^(O.l) ~ 0.3082, while for 
v = 0.01, this gives afit(O.Ol) ~ 0.0325 to be compared 
with, say, a^(O.Ol) ~ 0.0341. Note also that the an- 
alytically expected value in the — > limit is given 
by a = Vl2 v + 0(y^) = 3.4641 + 0{v'^), where 
the analytical value \fVl for the a\ coefficient derives 
from the well known specific angular momentum of a 
test particle at the LSO. The slight difference between 
^analytical ^ ^ ^ 3.46410 and af = 3.27690 is proba- 

^ \Vhile writing up our results for publication we became aware of 
the recent work [ssll where similar fits are advocated (for similar 
reasons to those discussed here). Note, however, that ^381] does 
not focus as we do on the losses during the ringdown phase, 



bly due to the deviations from analyticity in i/, as — > 0, 
implied by the appearance of strange fractional powers of 
V (integer powers of v^^^) during the transition between 
the LSO and the plunge, see Ref. [H. 

Let us finally mention that, though we focussed here 
on nonspinning binaries, we intend to study, within the 
EOB approach, the dimensionless spin parameter of spin- 
ning black-hole binaries. Indeed, not only is the function 
a(rni, m2. Si, S2) a useful diagnostics for comparing ana- 
lytical and numerical results, but it has also an important 
physical meaning. If the cosmic censorship conjecture is 
correct, this function should always stay smaller than 1, 
even if the individual spins take their maximum Kerr 
values Si = ml, 82 = "^l- Both leading-order EOB ana- 
lytical results [7, '2a|, and recent full-scale numerical [1§| 
results have indicated that this is indeed the case. It 
would be, however, quite interesting to improve the EOB 
estimates beyond the leading order, and to compare them 
in detail with numerical results. 
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but during the entire plunge -|- merger phase (post-LSO). As 
explained above, it is our use of the analytical EOB description 
which allows us to define precisely the "ringdown phase" . 
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